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Abstract 



We investigate the initial value problem for the Einstein-Euler equations of general relativity 
■ under the assumption of Gowdy symmetry on T^, and we construct matter spacetimes with 

low regularity. These spacetimes admit, both, impulsive gravitational waves in the metric (for 
instance, Dirac mass curvature singularities propagating at light speed) and shock waves in the 
' ' fluid (i.e., discontinuities propagating at about the soimd speed). Given an initial data set, we 

establish the existence of a future development and we provide a global foliation in terms of a 
I . globally and geometrically defined time-function, closely related to the area of the orbits of the 

^jq' symmetry group. The main difficulty lies in the low regularity assumed on the initial data set 

which requires a distributional formulation of the Einstein-Euler equations. 

1^ ! 1 Introduction 

<^ 

■<sj- . Background 
O 



We consider matter spacetimes with Gowdy symmetry which, by definition, admit a two- 
parameter group of isometrics generated by two orthogonally transitive, commuting Killing 
' fields. (See Section |2l below.) Under this symmetry assumption, the initial value problem for 

the Einstein equations has received much attention in recent years, both in the vacuum case and 
in the matter case when the matter is governed by the Vlasov equation of the kinetic theory of 
gases. In the present paper, we are interested in the evolution of perfect fluids, and in the context 
of Gowdy symmetry, we aim at constructing a globally foliated, future development of a given 
initial data set. The main novelty in this work lies in the presence of shock waves which appear 
in the fluid and, in turn, generate curvature discontinuities propagating in the spacetime. 

Recall that a spacetime is a (3 + 1) -dimensional differential manifold M that is endowed with 
a Lorentzian metric g with signature (-, +, +, +), satisfying the Einstein field equations 

Gap^xTap, (1.1) 

where k > is a constant and all Greek indices lie in the range 0, . . . , 3. Here, Gap ■- Rap-{Ri/ /2)gap 
denotes Einstein's curvature tensor. Rap the Ricci curvature, and Ry^ the scalar curvature of the 
manifold. The stress-energy tensor Tap appearing in the right-hand side of l|l.lt describes the 
matter content of the spacetime which, for perfect fluids, reads 

Tap := ijl + p) UaUp +pgap, (1-2) 
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where the scalar field represents the mass-energy density of the fluid and n its velocity vector. 
The spacetime is assumed to be time-oriented and u is normalized to be a future-oriented, unit 
timelike vector field. We also assume the linear density-pressure law 

p--k^fi, (1.3) 

in which the constant k e (0, 1) represents the sound speed in the fluid and does not exceed the 
light speed normalized to be 1. 

Under the assumption of Gowdy symmetry made in the present paper and after introducing 
areal or conformal coordinates (see Section|3l below), the Einstein field equations take the form of 
a coupled system of nonlinear wave equations with differential constraints. Since the pioneering 
work by Gowdy O, (vacuum) Gowdy symmetric spacetimes have been extensively studied 
ll24[|ll[[T7l [7l[8ll5ll4l[T8l and Penrose's strong cosmic censorship conjecture [SS",!^ was eventually 
established by Ringstrom |3T1 [32 J . A generalization of these spacetimes that contain matter 
governed by the Vlasov equation was recently presented in |26ll2gl [Tll2l lT0ll33l . 

As far as compressible matter is concerned, the mathematical investigation of Gowdy-type 
spacetimes was initiated by LeFloch and Stewart [22J (see also [3J), who introduced a converging 
approximation scheme for the initial value problem and derived several a priori bounds in suitably 
chosen coordinates. Therein, it was found necessary to cope with weak solutions to the Einstein 
equations, understood in the distributional sense and containing propagating discontinuities 
(shock waves). In ||3l|22l, the authors established a local-in-time existence result in the class of 
spacetimes with bounded variation. The present work is a continuation of this work and is aimed 
at constructing a global foliation of such spacetimes. 

Recall that spacetimes with boimded variation were constructed by Christodoulou in his work 
16 J that settled positively the weak version of Penrose's cosmic censorship conjecture in the context 
of spherically symmetric spacetimes and for scalar fields. Recall also that Groah and Temple [15 J 
established a local-in-time existence result for spherically symmetric matter spacetimes. In such 
spacetimes, no gravitational waves are permitted and the matter equations are coupled with a 
differential equation accounting for (non-evolutive) geometrical features. 

The study of solutions of low regularity is motivated by the fact that these can arise from 
smooth initial data. In the case of the Euler equations without gravity, this is well known and it is 
physically plausible that adding gravity should not make a fimdamental difference. A proof that 
this loss of regularity occurs in general relativity was given in [30 1. That result concerned plane 
symmetric solutions of the Einstein-Euler equations, and they are a special case of the solutions 
studied in what follows. A direct comparison is made difficult by the fact the time coordinate 
used in Ii30l . a constant mean curvature coordinate, is different from that used in the present 
paper. 

Objectives of this paper 

The assumption of Gowdy symmetry made in the present paper allows us to address the prop- 
erties of dynamical gravitational degrees of freedom. We consider the initial value problem, and 
search for the spacetime as a future development of a prescribed initial data set. This amounts 
to prescribing an initial hypersurf ace endowed with a Riemannian metric and its second funda- 
mental form, together with the mass-energy density and current vector determined by the fluid 
on this hypersurface. Our main result provides global existence for the Einstein-Euler equations 
within a class of spacetimes with low regularity which may contain, both, impulsive gravitational 
waves and shock waves: 

• When a shock wave arises in the fluid, the fluid variables become discontinuous and, as a 
consequence of Einstein's equations, the spacetime curvature becomes discontinuous. 

• In addition, our theory allows for distributional curvature singularities propagating at light 
speed. Such waves are referred to as impulsive gravitational waves. 
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• Our framework also encompasses the situation where the fluid contains low density regions, 
even vacuum regions. 

The theorem established in this paper can be regarded as a statement on the nonlinear stability 
of (vacuum) spacetimes with Gowdy symmetry on when compressible matter is included. 
Indeed, our matter spacetimes can be made to be arbitrarily close (in a well-defined functional 
norm) to vacuum spacetimes, by choosing the total amount of matter energy on the initial 
hypersurf ace to be arbitrarily small. 

More precisely, by combining geometrical and analytical arguments, we establish the existence 
of a global foliation based on a geometrically defined time-function coinciding with the area of 
the orbits of the symmetry group. Our approach is motivated by pioneering works by Moncrief 
| |24| (Gowdy spacetimes) and Berger, Chrusciel, Isenberg, and Moncrief H) (vacuum spacetimes 
with symmetry) which established the existence of such a foliation for vacuum spacetimes. 
However, we bring in a conceptually new and mathematically challenging aspect in that we 
consider solutions that have very low regularity. In the class under consideration, many (high- 
order) estimates derived in |4J no longer hold and must be bj^assed. The regularity of our 
spacetimes is considerably lower than the one constructed in earlier works and, for this reason, 
our analysis encompasses a larger class of spacetimes. 

Concerning the regularity of the spacetimes constructed here, the following features should 
be stressed: 

• Natural estimates for the geometry. As we will show, it is natural to impose that the essential 
metric coefficients belong to the Sobolev space of functions which, by definition, are 
square-integrable, together with their first-order derivatives. This regularity is dictated by 
the energy-type functional associated with the problem under consideration. 

• Natural estimates for the fluid. As we will also show, it is natural to impose that the energy 
density and the current vector belong to the Lebesgue space of integrable functions, but 
no further regularity can be imposed on their derivatives. Our setting allows the fluid 
variables to contain vacuum states and the essential velocity component to approach the 
speed of light. 

For further results on the existence and qualitative properties of spacetimes with matter when 
very low regularity is assumed, we refer to LeFloch and Stewart |22ll23l . 

An outline of this paper follows. In Section|2l we present our assumptions and main results; 
see Theorem l2.1l In Section|3l we express the Einstein-Euler equations first in a general foliation 
and, then, in areal and in conformal coordinates. Next, in both the expanding and the contract- 
ing directions handled in Sections |4] and |5l respectively, we introduce suitable notions of weak 
solutions to the Einstein-Euler equations (Definitions 14.51 and 15. It adapted to the setting under 
consideration, and we derive analytic and geometric properties. Finally, having successfully de- 
termined the natural weak regularity assumptions required on the initial data and enjoyed by the 
solutions to the Einstein-Euler system, we conclude with the compactness framework established 
in LeFloch |20I and presented in Section[6]below. 

2 Global foliations of matter spacetimes 

Symmetry and regularity assumptions 

First, we need to specify the symmetry and regularity conditions characterizing the class of 
spacetimes (M, g) iinder consideration. Throughout we consider matter spacetimes satisfying the 
Einstein-Euler equations (|l.ll l- lll.3l l. We assume that they admit an Abelian isometry group 
with spacelike orbits generated by two linearly independent, commuting, spacelike Killing fields 
K, L whose twist constants vanish, that is, 

:= eapy6K"Lt^y^'K' = 0, := ea^ys^K'^hNyi' = 0, (2.1) 
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where ea^yb denotes the totally anti-symmetric elementary tensor. This latter condition expresses 
the assumption that the frame of covectors (^g{K, •), g(L, - fj is integrable in the sense of Frobenius, 
and was identified by Chrusciel [7J in order to single out (vacuum) Gowdy spacetimes [14J within 
the larger class of symmetric spacetimes. Note in passing that the scalars c^, must be constant 
throughout any symmetric vacuum spacetime, as was observed by Geroch |T2llT3l . Matter 
spacetimes satisfying the conditions above are said to have Gowdy symmetry. The topology of 
the manifold, denoted by M, must also be specified, and we assume here that the spatial slices 
have the topology. All of the above conditions determine the class of spacetimes with Gowdy 
symmetry on T^, under consideration in the present paper. In the vacuum case, these are precisely 
the so-called Gowdy spacetimes first studied in [14]. 

We are interested in the future Cauchy development of a given initial data set, and the 
global topology will turn out to be M = [cq, ci)xT^ for some cq < ci < oo, and the spacetime to be 
foliated by spacelike hypersurfaces along which a time-function (denoted below by t or t) remains 
constant. The main unknowns of the theory are the Lorentzian metric g describing the geometry, 
together with the scalar field fi and the vector field u characterizing the matter content (via l|1.2b ). 
In local coordinates x = {t, xf) {a = 1,2, 3), Einstein's field equations dLljl will be decomposed into 
evolution and constraint equations for the metric coefficients which should also be coupled with 
the Euler equations for the evolution of the fluid variables (see below). It will be convenient 
to express the Euler equations as evolution equations for the mass-energy density p and the 
momentum measured by an observer moving orthogonally to the hypersurfaces determined 
by the time-function. By construction, j is tangent to the leaves of the foliation, and prescribing 
the fields p, ; is equivalent to prescribing the fields /,/, u (see (|3.2| |, below). 

We now discuss the regularity of the spacetimes under consideration. Our regularity assump- 
tions must allow for propagating discontinuities and, in view of the expression of the energy 
functional associated with the Einstein-Euler equations (see Lemma [4. 11 below), it is natural to 
require that on each slice labeled by the parameter f, both fields p{t, ■), i{t, •) belong to the space 
0{T^). In short, we write 

p,yeL-([co,ci),Li(r3)). (2.2) 

All the Lebesgue and Sobolev spaces under consideration are endowed with the volume form 
induced by the Lorentzian metric. By definition, for almost every f the function p(t) is well- 
defined and integrable, and the function t i-^ llp(f/ OIIl'ct-') (i-e- the integral of the function p(f)) 
is boimded iiniformly for almost all t in any compact subinterval of [cg, Ci). As we show later, 
the regularity Il2.2b is sufficient to formulate the Euler equations in the sense of distributions. In 
addition, following the general theory of nonlinear hyperbolic systems [9. 19], we impose that 
p, /' satisfy certain entropy inequalities (see | |4.12|| , below) which select the physically meaningful 
solutions to the Euler equations. 

On the other hand, again in view of the expression of the energy functional we impose that 
on each slice the essential metric coefficients belong to the Sobolev space . In short, we write 

geLZ,{lco,c,),H\T^)). (2.3) 

It should be noted that our precise definition, given in Sections |4] and |5]below, are (slightly) more 
general and allow certain components of the metric to be less regular than . 

Due to the assumed symmetry, it will turn out that l|2.3|l implies that the metric coefficients 
remain locally boimded, that is, g e L^^, and our regularity condition suffices to express Einstein's 
equations in a weak sense. Indeed, the Riemann curvature of g and all of its traces and, therefore, 
the Einstein tensor are then well-defined as distributions [21J. Note finally that the regularity 
described here may in principle depend on the foliation under consideration, and may not be as 
geometric as one may wish. 

We refer to such a set (M, g, p, j) having the regularity | |2.2|I -| |2.3| I as a finite energy spacetime 
with Gowdy symmetry on T^. 
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The initial value problem 

To formulate the initial value problem we prescribe a Riemannian metric ^ on together with a 
2-covariant symmetric tensor field k, both of them satisfying the assumption of Gowdy symmetry 
on T^. By definition, the Lie group acts as an isometry group on the torus generated by two 
(linearly independent, commuting) vector fields, 

'Cxg = 'CYg = o, [x,y] = o, 

satisfying the condition of vanishing twist constants l|2.1) l. We also prescribe a scalar field p and 
a vector field defined on T^, also satisfying the Gowdy symmetry assumption, and we require 
the regularity 

geH\T% keL^{T% p~jeL\T\ (2.4) 

(As already pointed out, we actually cover slightly more general data in which certain components 
of the fields ^, A: are less regular.) Finally, we assume that these data satisfy Einstein's constraint 
equations 

R + (tr l<f -f\= Ik p, tr (v fc) - v(tr k) = k], (2.5) 

where V and R denote the covariant derivative operator and scalar curvature of (T^, 'g), respec- 
tively, and the trace and norm are determined from the metric "g. Under the regularity l|2.4l l, the 
curvature terms R, tr (Vfc),and V(trfc) are well-defined as distributions [21J, while the remaining 

terms in | |2.5l l belong to L^(r'). Of course, | |2.5l l actually implies that the scalar field R is more 
regular and, actually, integrable in space. In fact, we will not use directly the definition from 
IJlJ; instead, later in this text, a complete definition of weak solutions appropriate to the problem 
under consideration will be introduced. 

We refer to {g,k,^,fj as a finite energy, Gowdy synunetric, initial data set on for the 
Einstein-Euler equations. 

We are now in a position to state the main result of the present paper. Recall first that an 
important quantity associated with Gowdy-symmetric spacetimes is the area function R (defined 
up to a multiplicative constant) of the orbits of the symmetry group generated by the Killing 
fields. Since the metric coefficients are only in H^, the gradient VR may only be defined almost 
everjrwhere for the Lebesgue measure. However, using the Einstein equations we will prove 
later (cf. Proposition l5.2i below) that VR is actually continuous. In addition, a standard argument 
showing that VR is a timelike vector field remains valid under our regularity condition and 
even for compressible matter. This condition on the area function can be expressed explicitly in 
terms of the initial data: when the spacetime metric is sought for in conformal coordinates (see 
(|3.261 , below), the initial data prescribed for the time derivative of the area function (denoted by 
Rq below in Section |5]l should be everywhere positive (expanding case) or everywhere negative 
(contracting case). 

This property allows us to distinguish between two cases, depending whether the spacetime 
is expanding or contracting, that is, the area of the orbits of symmetry is increasing or decreasing 
toward the future, respectively. Furthermore, without genuine loss of generality we assume 
that the initial hypersurface has constant area function, which imposes certain restrictions on the 
initial data set but, experience with (vacuum and regular, at least) Gowdy symmetric spacetimes 
tells us that this is not a restriction on the class of spacetimes themselves, only a minor restriction 
on the choice of the initial hypersurface. 

The following main theorem encompasses both the expanding and the contracting cases. 

Theorem 2.1 (Einstein-Euler spacetimes with Gowdy symmetry on T^). Let (g,fc,p, be a finite 
energy, Gowdy symmetric, initial data set on for the Einstein-Euler equations, and assume that these 
initial data have constant area, and are everywhere expanding or everywhere contracting. Then, there 
exists a finite energy, Gowdy symmetric spacetime (M, g, p, j) satisfijing the Einstein-Euler equations 
(|l.l)l - l|1.3|l in the distributional sense, and the following properties hold. The manifold {M,g,p,j) is (up 
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to dijfeomorphisms) a Gozvdy-symmetric future development of(^g,k,p,j^, which is globally covered by a 
single chart of coordinates t and {6,x, y) e T^, with 



M={ 



0) / < Co < t < cx)| X T^, expanding case, 
^{t,d) I cq < t < ci < 0| X T^, contracting case. 



Here, ci e (cq, 0] is a constant, and the time variable is chosen to coincide with the area of the surface of 
symmetry in the expanding case, and with minus this area in the contracting case. 

The above theorem relies on a notion of weak solution which will be described in full detail in 
this text (in Sections 4 and 5); importantly, this definition incorporates the initial data in a weak 
sense as well. 

It is interesting to ask if ci can be taken to be equal to zero. Note for comparison that if the 
perfect fluid occurring in the theorem is replaced by collisionless matter described by the Vlasov 
equation then the analogous theorem does hold with ci = 0, as was established in |35| . For a 
fluid with a linear equation of state there are solutions for which the theorem does not hold with 
Ci = as follows from the discussion in Section 12.2 of |29]. At the same time it seems, as will be 
discussed in more detail in Section|3l that these solutions are exceptional in this regard. 

The above statement provides a global foliation of the constructed spacetimes. Due to the time 
irreversibility of discontinuous solutions to the Euler equations, compressible matter spacetimes 
can only be defined in the future of the initial hypersurface. This is in contrast with vacuum 
spacetimes and Vlasov spacetimes which can be defined in both time directions. 

An important and very challenging open problem is to establish the strong cosmic censorship 
(that is, the inextendibility of the future Cauchy development) for the spacetimes constructed in 
Theorem l2.1l In the expanding case, the future inextendibility of the spacetime follows from the 
fact that the area function tends to infinity along every future timelike direction. The generic 
inextendibility of (vacuum) Gowdy spacetimes in the contracting direction was established in 
major contributions by Chrusciel, Isenberg, and Moncrief |8| (polarized case) and Ringstrom IISTl 
|32J, and recently also established by Dafermos and Rendall for the Vlasov model [10] (cf . Theorem 
4.1 therein). 

Observe that the uniqueness issue is not addressed here, since within the functional frame- 
work under consideration such a result is not even known for the Euler equations in the flat 
Minkowski geometry. The entropy inequalities are however fully motivated, and are known 
to imply uniqueness for certain restricted classes of initial value problems or smaller classes of 
regularity of solutions; we refer the reader to [9. 19 1 for details. Note also that although the areal 
foliation constructed in Theorem l2.1l is geometric in nature, it is tied to the Gowdy-symmetry as- 
sumption, and it would be interesting to investigate the existence of a global foliation by constant 
mean curvature (CMC) slices. 



3 Einstein-Euler equations 

A (3 -I- l)-decomposition of the Euler equations 

We present here a formulation of the Euler equations in a divergence form, which we will later 
use to define a notion of weak solutions. All calculations in the present section are carried out in 
the smooth setting. Recall that the Bianchi identities for the geometry implies the Euler equations 
for the fluid 

v"rf = 0, (3.1) 

in which the energy-momentum tensor determined by (|1.2|l . We will eventually express these 
equations in local coordinates adapted to the Gowdy symmetry but, as a preliminary step, we 
consider a general foliation by spacelike hypersurfaces, determined as the level sets of a time- 
function f. Choosing a vanishing shift vector and introducing the future-oriented unit normal 
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n", the lapse function N > 0, vanishing shift-vector, and the second fundamental form k"f^ to the 
foliation, we can write 

k{X,Y) := giVxXn) = -g(Vyn,x). 

Here, for each hypersurf ace, X, Y denote arbitrary vector fields tangent to any given hypersurface 
(with y arbitrarily extended to a neighborhood of that slice). 

Using the normal n we determine the projection operator h"l^ := g"^ + satisfying the 
obvious conditions h"^na = and h"t^Xa = X'^ whenever UaX" = 0. We decompose the matter 
tensor T"'^ into its normal and tangential components, as follows: 

p := T^^Hatip, f := -Tl^^h'-^ny, S"^ := V^^hl, (3.2) 

so that T"^ = pn"n^ + fn^ + fn" + S"^*. The scalar p is the mass-energy density measured by an 
observer moving orthogonally to the slices of the foliation and / is its momentum vector which 
is tangent to the hypersurf aces. Using the expression l| 1.21 1 of the energy-momentum tensor, we 
find ^ 

p = (p + p)[uan") -p, f = -ip + p){uiinf^)u^'h"y. 

Note that T^hiatip = p, so that p is the mass-energy density measured by a (Lagrangian) observer 
moving with the fluid. 

In coordinates (t, x") {a = 1, 2, 3) adapted to the foliation so that f is constant on each slice, we 
have (gafi) = -N^ df + (gab) and 

1 

no = -N, n„=Q, kab = -^dfgab, 

p = N^T°°, f = 0, f=NT°\ (3-3) 

5°" = 0, = 0, S"^ = r%h^^. 

Since ifria = tfifio = -Nu° and -1 = u"Ua = -N^(u'')^ + gabu'^u^ , we thus find 

p = {p + p)N^{u°f - p, f = (p + p)Nu°u\ 
S"'' = {p + p) u''u''+p^\ 

Now, writing the Euler equations l l3.Hl as 

= V«Tf = [p (n„ + Uy + fn„ + 

= nPV" (pria + U) + V"si + (V«/)«a + ^"n^ {pn^ + ]„) + f^n^, 
and decomposing them into normal and tangential components, we find 

= -V" (pn„ + ;„) + npV"si + naUpiV^f), 

= h'yV'jyna + h^S^Sl + h^y"ny ipn„ + ja) + fV'na. 

At this juncture, we observe that in order to define a concept of solution to the Einstein-Euler 
equations as we will do later, it is important to keep these equations in divergence form. The 
suitable form of the equations consistent with the definition of distributional geometry presented 
in [21 1 will be obtained provided we take into account the volume element associated with the 
Riemannian metric ^^'^ induced on the slices, as well as the volume element associated with the 
Lorentzian metric g. With this aim in mind, we introduce the rescaled matter fields defined as 
follows: 

p:=cvp, a; := (det (3)^)1/2, 

~ - . , 3.5) 

j := N i := N CO j, S := S := co S. 



(3.4) 
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The scalars n and p are rescaled in the same way as p, with an obvious notation. Observe that 
scalar, vector, and 2-tensor fields are scaled differently. 

After some further calculations, the Euler equations (|3.H l take the form 

dtp + daf = El, 

<9,f +wP)Vb(w-iS''^) = E«, 

where ^^'V is the connection induced on the hypersurfaces, and 

El := N-H-daNj" +kaiS"% 

„ (3 7) 

E^ := N-^{dtNf + dbNS"'') + 2Nk"tf - pNVN. 

This completes the derivation of a general (3 + l)-decomposition of the Euler equations. Note 
that the equations I l3.6b have a divergence form. 

Now, under the Gowdy symmetry assumption the Euler equations simplify, provided we 
use coordinates x - {t,x^,x^) (with B = 2,3) adapted to this symmetry, so that the coordinates 
{x^,x^) e describe the 2-surfaces spanned by the Killing fields K,L. In agreement with the 
derivation of I|3.6t -l l3.7|l , our choice of coordinates {t,x^) is made so that goi ~ (see the metric 
expression 1 13.131 and | |3.26|I , below). Furthermore, imposing 

= 1/3 = (3.8) 

is compatible with our symmetry assumptions. So, recalling that u" is a future-oriented, time-like, 
unit vector, we find 

NV)'-^ii("')' = 1- (3.9) 

The velocity vector is determined by the single function u". 

The (non-vanishing) components of the matter tensor then take the form (with B, C = 2, 3) 



p = {^t + p)N\u°f - p, f = {ix + p)Nu°u\ 
S" = (fi + p) (u^f + p g", S^^ = 0, S^'= = 

From their definition, we easily obtain the Christoffel symbols 



and we are in a position to compute explicitly a key term in H3.6|l : 

a;(3)v,(a;-iS«') = ^iS" + TIS'^ + F^^S"^ - cv-'VycvS'" 
= 5iS" + Fj,S" + Fi^S^^. 
Recalling that [i and p are defined by rescaling as p in | |3.5| |, we arrive at the formulation 

dfp + dif = El, 

(3.10) 

dtf + (9iS" = E2, 

the right-hand terms being given by 

El :=JV (fcii Qi + p) {u^f + (tr fc) p) - N diN (p +p} u°u\ 

E2 = Ll ~N (dtN{p+p)uV +diN((p+p)iu^f +pg^^)) 

^2 (3.11) 

+ 2 NV^ Qi + p)u^u^-p diN. 
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Denoting by ^^'g the determinant of the surface of symmetry and observing that 

^"^■^- = 5^piiy*<*'"'^)- ^^^^^ 

we conclude that in order to express the right-hand sides Ei, E2 of the Euler equations | |3.10|| , we 
need compute only the zero-order terms 

N, gn, det(3)g, det^^)^, 

and the first-order terms 

dtN, diN, dign, 5i(detP)g), diidet'^''^ g), k\. 
Based on ll3l0l l- (|3l2l l, we are now in a position to write the field equations in local coordinates. 

Areal coordinates 

Recall that this discussion is carried out in the smooth class; weak solutions in areal coordinates 
will be discussed in Section 4, below. We begin with the so-called areal coordinates {t, 6, x^, x^) 
in which the metric is determined by four scalar functions U, A, j], a, depending on the variables 
(f,0)only: 

g = e^*''-") {-a df + dd^) + e^"(dx^ + A dx^f + e'^^ f {dx^f (3.13) 

and the variables d, x^, x^ range over the interval [0, 1ti\ and are 27i-periodic. The coordinates t, 6 
parameterize the quotient manifold M/T^, while d/dx^ and d/dx^ denote the Killing fields and 
x^,x^ are coordinates on the torus T^. By construction, the area of the two-dimensional spacelike 
orbits of symmetry coincides with the time variable, since 

'^""^{Ae^" e-2"t2+A2e2"J = ^ • 

The fields of 1 -forms 

eo ■- a^^^e*''-"' dt, ej ■- e'^'^dd, := {dx^ +Adx^), £3 ■- e'^tdx^ 
determine an orthonormal frame, whose associated dual frame is 

-i/2-(,,-u)^ e-(')-(.r)A p-'^— e^'-( -A — + —\ 

dt' de' dx^' dx^ dx^r 

Returning to the matter variables and recalling the "projected" energy-momentum tensor S"^ 
introduced in (|3.2|l , we extract its essential components 

Pi -S(ei,ei) = e^''-2"S\ 

Pi ■■= 5(62, 62) = e^" (t22 +2AT^^+A^ T^^), (3.14) 

P3 := 5(63, £3) = f P23 := S(e2, £3) = tT^^ + tA 

From now on, we use subscript indices to denote partial derivatives with respect to t, 6. After 
a tedious calculation, Einstein's evolution equations extracted from l(Ll| take the form of three 
nonlinear wave equations for the metric coefficients U,A, rj 

Utt-aUee-(-] + ^)Ut + ^Ug + Q"(a,0 + aU'-'it), 

An - aAgg ^ + ^)At + ^ Ag + Q^(a) + atn^, (3.15) 

(Xf; CXq 1 ^'q 

ntt - ariee = i^rit + — rfe + ij^agg - — + □''(«, t) + a lV\t), 
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in which the lower-order metric terms are 

Q"(a, t) := — (Af - a A^), Q^a) := -4(U,A, - a UeAe), 

Ml 

GL\a,t) := i-Uf + a + ^(A? - aA^), 
and the lower-order matter terms are 

:= (p - Pi + P2 _ P3), := 2 e^O-zu) p^^^ 

W{t) := -iA2e2(')+i^)p2 - e2('i-")p3 - jAe^'^P23. 

On the other hand, Einstein's constraint equations extracted from provide us with three 
first-order differential equations for the metric unknowns r] and a, that is. 



-T], = Uj + aUl + —^ (Af + aAl) + e^^'^-^'^a p, 



1 e^" 
yT]f = Ul + aUl + 

yrie = - ^J^^e + 2 LJ^Ue + ^ A,Ae - e3W-")ai/2yi^ (3.16) 



yat = -2a2e20,-u)(p_p^)_ 
To finally express the Euler equations, we need the zero-order metric quantities: 

as well as first-order ones: 

tr fc=-a-i/2e-(''-"'(i+7],-!Jf), 

fci = giifcii = _a-i/2e-(')-y)(;jj _ tJ,). 
These expressions allow us to rewrite the equations (|3.111 in areal coordinates: 

(p}t + {J)e =U, 
+ (S")e =E2, 
with now 

El = - e^^'^-''\r]t - Ut){Ti + v){u^f -(\ + nt- Ut)p 
-lde(ae'^'^-^yji + p)u'^u' 

U=\ d{a e^C'-")) CI + p) m°m1 + i o)e(a e^Ci-u)) ((^ + p) (^1)2 + pe-2(']-u)) 

- a (r]0 - Ue) (e^CJ-^^) + p)(m1)2 + p) 

- 2ae2('?-u)(,^^ _ LIf) p) uOi/i - ip(ae + 2a(r]0 - iJe)). 

Finally, to "close" the system it is necessary to have explicit expressions of the conservative 
and flux variables 'p, f, S^^ (arising in | |3.17| |) in terms of the energy density and the velocity 
function 



and 



(3.17) 



V := a-y^ 
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Noting that a {u°f - {u^f = e-^Oj-Lf), we obtain 



g2(,,-U) („0)2 ^ g2(,,-U) (yl-)2 ^ 



a(i-y2)' ' ' i-y2- 

In view of p = /.i + (p + i.i)e^^'i~^ (u^)^ , the rescaled mass-energy density reads 

p = fe(''-^)(^ + (p + p)^^). (3.18) 

For the momentum we find = e"^''"^^ (/.i +p)y(l - so that the rescaled momentum variable 
reads 

7 = tai/V'-"(^i + p)^^. (3.19) 
Finally, from S^^ = tae^^'^~^^{i.i + p)(u^)^ + pg^^^ we obtain 

S" = te(''-")a ((H + p)^^ + p). (3.20) 
Consequently, we can express Ej and E2 in terms of fi, V and obtain 

^i=- e''-" (f ('?f - Ut){^ + p) + (1 + t(ij, - !Jf))pJ 



and 



1:2 =^ a-i/^e-J-Lf (a, + 2a{r]t - U,)) (p + p) ^— ^ 

+ 2 («9 + ^"(^e - !i9))((p + p) J— ^ + p) 
- 2f a (T]e - Ue) e''"" ((p + P) + p) 

-2fai/V'-"(7],-LJ,)(p + p)^^, 
which completes the derivation of the Euler equations in areal coordinates 

(te("-^"(p + (p + p)^)) +/te"-"ai/2(^ + p)-Ji-^j .El, 
|fe"-V/2(p + p)^Jl^J +|te("-")a((p+p)^^+p)) 



(3.21) 



Finally, assuming the linear equation of state ( 11.31 1 we conclude that the fluid variables p, V 
satisfy the nonlinear hyperbolic system of two balance laws 



(3.22) 



|^,.-V/2 ill^p)^ + (^.-r^^a ^p|^ =f p."-"2:^', 
with coefficients depending upon the metric functions U, i], a, and 



(1 + k^)V 
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and 

^2 =«"'^' (y - <^{nt - ^^))-Yz4^ + ( Y - <ne - Ue)) YT^- 

As we will see later the second equation in (|3.22l l should still be multiplied by a~^, leading to the 
final form of the Euler system 

+ , .1/2(1 + Si 

(3.23) 
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te^ a ' —u\ +ue'"" u =te'"" — = — u 
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with 



and 



-01 := ^ 



+ {r]t - Ut) Q^ + V^) + a^l^ + 'Je " ^e) (1 + k^W 



-S2 ■.=a-'l^ + 1' - + l)V + (^+^g-Ue){V^ + k^). 

Observe that the metric coefficient A does not arise in the Euler equations. 

It remains to compute now the expressions of Pi, P2/f 3/^*23 required in the right-hand side 
of | |3.15|I . Computing the coefficients of the inverse of the metric, specifically the coefficients 

g22 ^ r'^e^^A^ + £-2", ^23 ^ _f-2g2u^^ ^33 ^ r'^e^^ , we find 

^l^T^^T^^'' ^2 = ^3=? = ^"^/, P23 = 

and, therefore, by recalling that p = (fi + pV'^)l{l - V"^), 



n''(f):=-e'''fc'(e-2" + 4e'")fi. 



(3.24) 



t2 

Note also that the expression of the right-hand of the equation for at 

y(l/a)f = 2(l-fc2)e2(TU)^j 
(l-F)(l-y^) 



(3.25) 



Conformal coordinates 

Recall that this discussion is carried out in the smooth class; weak solutions in conformal coordi- 
nates will be discussed in Section 5, below. One may also use coordinates (x, 6, x, y) in which the 
metric has the "isothermal" form 

g = e^f'i-") (-dx^ + dd^) + e^"{dx + A dyf + e'^" dy^ (3.26) 

and thus depends on the four scalar functions rj, U, A, and R of the variables (x, 6), only. The area 
R of the surfaces of symmetry is now an independent unknown function. Formally, the function 
a is replaced by 1, and the variable t is replaced by R, and one therefore still has to determine four 
metric coefficients. 
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In comparison with the areal coordinates, the Einstein-Euler equations take a somewhat 
simpler form in conformal coordinates. There are now four evolution equations: 



where 



U„ - Uee = + ff", 
A^-Agg =Q^(1) + R^^ 
r]„ - r]ee = + ff'', 
R„ - Ree = ff"". 



:= ^(RtAx - RgAg) + Q^(l), Q'' := Q''(1,R), 



(3.27) 



and 



■- n"{R), := n^, n'' := W{R), 

Hence, the lower-order terms are given by (essentially) the same expressions as in Section |3l the 
i-derivatives being replaced by T-derivatives, and a, t replaced by 1, R, respectively. In conformal 
coordinates, there are only two constraint equations 

^ = ^iVrRr + noRe) - iU', + Ul)-y (A? + Al) - p, 

R 1 (3.28) 

-f- = ^(VrRe + VeRr) - 2 UMe - ^ A.Ag + e^Oi-t-Oyi, 

which can be regarded as first-order differential equations for Rg. Furthermore, the Euler equa- 
tions I I3.23I I can similarly be written in conformal coordinates. 



Spatially homogeneous solutions 

In this section the equations will be expressed in areal coordinates. A special type of solutions are 
those for which all variables are independent of the coordinate 6 and so all the evolution equations 
reduce to ordinary differential equations. These are what are known as spatially homogeneous 
solutions. They admit an action of a three-dimensional Abelian group by symmetries and so in the 
usual terminology of general relativity they are of Bianchi type I. It follows from the momentum 
constraint that - 0. By a linear transformation of the spatial coordinates the initial data A(0) 
and (9fA(0) for the variable A can be set to zero. It then follows from the evolution equation for 
this quantity that A vanishes at all times. The transformation is naturally defined on the universal 
covering of M and need not project to M itself. This fact is, however, irrelevant for the study of 
the dynamics. Thus in the homogeneous case it may be assumed without loss of generality that 
A vanishes identically so that the metric is diagonal. The Bianchi type I solutions of the Einstein 
equations can be determined explicitly provided a new time coordinate T is introduced which is 
a function of the areal time coordinate f . The explicit form of the metric is (see for instance [3ll, 
p. 199) 

- a^'dT^ + T^P^a^^de^ + T^P'a^'^^dx^ + T^P'a^i'dy^, (3.29) 

where the p, and qi are constants for / = 1, 2, 3 which satisfy the Kasner relations pi+pi+Ps = 1/ 
Pj + P2 + P3 - 1/ ai^d - \ - Pi- The function a{T) is defined by a^"*^^ = a + nP-T^~^ where a and 
m are positive constants. Comparing this with the general form of the metric in areal coordinates 
reveals that 
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This can not be solved explicitly so as to get T as a function of t and thus an explicit expression for 
(|3.29t in areal coordinates. It is nevertheless clear that the mapping from T to t can be inverted 
uniquely. The initial singularity occurs at T - 0. The parameter pi belongs to the interval [—1/3,1]. 
Provided pi < 1 it can seen that t goes to zero as T ^ and this corresponds to the situation ci - 
in the theorem. 

An exception occurs when pi = 1 since in that case Ci is strictly positive. The behavior of the 
geometry as f ci in this case has been analysed in l29l . The Kretschmann scalar RajSydR"^'^"'^ 
tends uniformly to infinity so that this is really a curvature singularity. On the other hand the 
metric can be extended to t = Cj in such a way that it remains continuous and non-degenerate. 
Moreover, with repect to this extended metric the hypersurface t = cj is null. This is what is 
known as a weak null singularity. Note that the above calculations show that within the class of 
spatially homogeneous solutions the exceptional solutions are nowhere dense. It is not clear what 
will happen with more general inhomogeneous solutions but all available information indicates 
that further solutions of this kind will be at best rare and possibly not exist at all. 



4 Expanding spacetimes 

A weak formulation of the Einstein equations 

We now discuss weak solutions to the Einstein-Euler equations, by relying on areal coordinates 
and considering the case of expanding spacetimes, in which the area function R is increasing 
toward the future from an initial value denoted by Cq. Hence, the initial value problem is posed 
with data prescribed on the hypersurface f := Cq > 0. Instead of /], it will be convenient to use the 
new variable 

1, 

V := 7] + - ma, 
/ 2 

so that the geometry is now determined by the four scalar functions U,A,a,v. The natural 
regularity imposed on these functions is determined in the following lemma. 

Lemma 4.1 (Energy estimate). In areal coordinates the total energy 

E = E{t):= f Za'^l^de, 

£ := (Ul + aUl) + ^ e*" (Af + a A^) + e2(v-u) 
is a monotone decreasing function of the variable t > cq and, more precisely, satisfies 

4-£ = -- r 1)de, 

dt t Js, 

© := ^2 + ^ e4U^l/2^2 + 1 ^-l/2g2(v-U)(p + p^) > 0^ 

The underlying structure is analogous to two quasilinear wave equations for the variables 
U, A, taken to have finite total energy 

£ {{a-^'^Uf + uD + 1^ e^u (^-1/2^2 + «i/2^2))^e^ 

coupled to the "first" Euler equation for p taken to have finite total mass-energy 

r pe2(v-u)^-i/2^0_ 

Even in the class of weak solutions, the relations in the above lemma hold as equalities in the 
sense of distributions. 
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In view of the expression of £, as long as a, a ^,U,v remain bounded, i.e. 

\a\ + \a\-^ + \U\ + \v\ < 1, (4.1) 

and provided the initial energy E(co) is finite, which we now assume on every compact subset of 
[cq, Lemma 141] allows us to control the norm of the functions Ut, Ug,At,Ag on every slice 
of the foliation, as well as the norm of the mass density, from the same quantities evaluated on 
the initial slice. Recalling the expression (|3.18l l of the density in terms of the fluid variables V, 
we see that the energy estimate imposes 



Ut, Ue,At,Ae e L^{Ll), M := e L^{Ll), (4.2) 



at least on every compact set in time. In view of ( I3.15K we see that U and A satisfy (second-order) 
wave equations and that two initial conditions must be prescribed for each variable, that is, 

Utico) = UeLl, U{cq) = UeHl, 

Atico) = Ae Ll, A{co) = A e H^. 

The initial data for v is then computed from the constraint equation l|3.16|l ; see l|4.13l l below. 
Next, in view of the constraint equation (see l|3.16l l) 

Ive = 2 UtUo + ^e^^'AtAo - e^^^-^^a-'l^f, (4.3) 

and since |/^| < p we deduce that vg £ ig- Similarly, for the time derivative we have 

1 p'^^ 

-vt = U^ + aUl + ^(Af + aAl) + e^^^'-^'^Pi > 0, (4.4) 

and using < p e T^^C^t) we obtain Vf e L^(Ef). It is natural to assume the regularity 

vt,vgeLr{Ll). (4.5) 

Instead of these differential constraints, we can also use the evolution equation for v which follows 
from | |3.15l l, in which case two initial conditions are required on the variable v 

v,(co) = ve4, v(co) = veWg'\ 
Finally, the constraint equation on the function a (see (|3.161 again) shows that 

e Lnil), (4.6) 

1 — V- 

n can be imposed ( 

a on the initial slice, that is. 



1- V2 

while no regularity condition can be imposed on the spatial derivative ag. One should prescribe 



a(co) = a e Lg. 

Having identified the basic functional spaces of interest for each variable, we now consider 
the algebraic structure of the equations. We must reformulate the Einstein equations in a form 
that makes sense under the limited regularity above, only. The following key observation can be 
checked by direct calculation from | |3.15l l. 

Proposition 4.2 (Weak form of the Einstein equations). The evolution equations for the functions U, 
A, and v take the form 

( a-y^At) - fr^ a^^ Ag) = -4ri a-y^ (UfA, - a UgAg) + a^^ n^, 
{t a-y\v + (1/2) log a) - {t a^^ vg)g = t a^^ {-Uf + a Ul) 

;iA^-aAl) + a''^n^, 



4fai/2 
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ivhere W - TP. 

Finally, we supplement H.Ji with an equation for the function a 

a{t, 6) = a id) exp ( - 2(1 - A:^) ^' (e2(v-(j)^ ^-^ _ v^)yt', 6) dt'), (4.8) 

where a is a prescribed data in satisfying Zv > 0. Clearly, since |Pi| < p the function a is 
globally bounded: < a < a . 

Importantly, in the proposed setting the Einstein equations I l4.7b now make sense under the 
regularity conditions (|4.2b -| |4l6l l. For instance, a~^^'^Ut is the product of an U° function by an 
L^(Lg) function, and its derivative is defined in the distributional sense. 

Remark 4.3. 2. The function v may he also determined by the constraint equation, that is, 

v{t,e)= f (2tUtUg + ^e^^^AtAg-t{l+k^)e''-^MV]{t,e')de', (4.9) 
Jo 2r 

which is an integral equation for the unknown v. 

1. We emphasize that, instead of the equation for v, the original evolution equation for rj can not be used 
to derive a priori estimates, since it would require second-order derivatives of a - which can not be obtained 
in general for the non-vacuum Einstein equations. This is clear from the equation in at: it involves the 
fluid density which generically contains jump discontinuities and, therefore, is not sufficiently regidar to 
allow us to control second-order derivatives of a. 

Finite energy solutions 

Let us further discuss the matter variables and introduce a suitable notion of solutions to the Euler 
equations. Recalling that M - p/{l — V^), and using the variable v instead of /], the equations 
(|3.231 become 



(4.10) 



t a-^^^e"-"{l + k^ V^)Mj + \^t e""" (1 +k^)VMj =t e^'^a-^'^ Si M, 
{te"-^ a-^{\+k^)VlsA^ +[te'-^ a'^^^ik^ + V^)M^ = t e"-'^ a'^'^ S2 M, 
with 

t 2a 

and 

-S2 :=a-^^2 (vf - Ut){k^ + l)V+{vg- Ug) {V^ + k^)- 

The selection of physically admissible solutions to the Euler equations will be based on the 
conservation law of the number of particles 

VaN" = 

satisfied by sufficiently smooth solutions, where the vector field N" has the form N" = ru" with 

dr _ dp _ k^ dp 
r ' p + p{p) k^ + \ p ' 

hence 



16 



For simplicity, we have assumed here that a single thermod5mamical variable suffices to determine 
the state of the fluid under consideration. According to this approximation, for the theory to 
encompass discontinuous solutions, we must relax the above equality and impose that solutions 
satisfy the inequality 

- VaN" < 0, (4.11) 

which we refer to as the fundamental entropy inequality. 

More generally, an infinite list of (mathematical) entropies is available. 

Definition 4.4. A vector field J is called an mathematical entropy flux for the Einstein-Euler equations 
if every smooth solution to these equations satisfies the balance law 

V„J« = S, 

where J" and 9 both contain lower-order terms only, that is, zero- and first-order terms in the metric 
variables and zero-order terms in thefiuid variables. It is said to be convex if the component 3^ is convex 
in the conservative variables of the Euler equations (that is, in the variables (p, j^)). 

Physically admissible solutions should then be characterized by the entropy inequalities 

V„ J" < g (4.12) 

understood in the distributional sense, for every convex entropy pair that vanishes on the vacuum. 
As presented in [20], the family of convex weak entropies can be described by an explicit formula 
which can be used to establish a compactness result for sequences of solutions. 
We are now in a position to state precisely our notion of solutions. 

Definition 4.5 (Finite energy solutions). 1. A finite energy initial data set for the Einstein-Euler 
equations is a family of measurable functions U, U,A,A,v,v,a, M, V satisfying 

inf a > 0, M > 0, \V\ < 1, 

SI 

together with the regularity conditions 

U,A,eL^{S^), U,AeH\S^), v,MeL\S^), 



VeW^'\S^), a,VeL°°{S\ 
and the constraint equation 

nd . _ _ 

v{0) = J [icoUUe + ^e'^^ AAe-coe~-^{l+l<?-)MV){e')de'. (4.13) 

2. A finite energy solution to the Einstein-Euler equations (in areal coordinates) is a family of 
measurable functions U,A,a,v,M, V defined on [cq, oo) x and satisfying the following conditions: 

• The functions a,M, V satisfy (for all ti > cq) 

inf a > 0, M > 0, |V| < 1. 

[co,ti]xSi 

• The geometric coefficients U,A, a, v satisfy the regularity conditions 

Ut,At,eL^J[co,+oo),L^(S% U,A e L^J[co,+oo),h\S% 
vt,ve e L^,([co,+oo),Li(Si)), MSi,MS2 e Ll^{[co,+^) x S^). 

• The functions U, A, v, a satisfy the evolution equations ||4.7|| and the constraint equation (|4.8|l in the 
distributional sense. 
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• The fluid variables M, V satisfy the entropy inequalities 



in the distributional sense for all convex weak entropy flux J" to the Euler equations. 
• The prescribed initial data are assumed in the sense of distributions. 

To be more explicit, let us for instance state the first equation in ll4.7t in the distributional sense: 

Jca JQ 



(ijjtta-^'^Ut-4>6ta^'^Ue + 4'-—^{A^-aAl) + ipta^'^n'^)dedt 

Jca Jo ^''^ 

ip{cor)coa~'^'^Ude = 




for every smooth, compactly supported and spatially 27T-periodic function i/^ : [cq, oo) x [0, In] —> 
K. Observe that the initial data is taken into account in this statement. 

Under the condition in the above definition, we will also use the phrase "spacetimes with finite 
energy". It is important to observe that, under the regularity assumptions under consideration 
all the terms arising in the equations under consideration do make sense, at least as distributions. 
In particular, the right-hand sides of the equations | |4.10| | belong to O on every compact subset. 



5 Contracting spacetimes 

A weak formulation of the Einstein equations 

We continue our discussion of weak solutions to the Einstein-Euler equations, using now confor- 
mal coordinates in the case that the spacetime is contracting. Recall that, in conf ormal coordinates, 
the metric is determined by the four functions U, A, i], R, while the fluid is determined by the func- 
tions M, V. 

It is convenient to pose the problem on a hypersurface of constant negative tq. 

Definition 5.1 (Finite energy solutions in conformal coordinates). Given tq < 0, an initial data set 
with finite energy /or the Einstein-Euler equations in conformal coordinates is a set of functions 



R e L'^iS^), R e W^''^(S^), inf R > 

SI 

and 

MeL\S^), yeL°°(Si), M > 0, |y| < 1. 

Given an initial data set as above, a set of functions U,A, r],R,M, V defined on some interval [tq, ti] c 
[to, 0) is called a finite energy solution to the Einstein-Euler equations in conformal coordinates if the 
following conditions hold: 

1. The functions have the regularity 

U„A, e L°°([to, Ti], L\S^)), U,Ae L^([to, ti], W'-^S^)), 

R, eL°°([To,Ti],L°°(Si)), ReL-([To,Ti],Wi'°°(Si)), inf R>0 

[To,Ti]XS1 

and 

M e L°°([to, Ti], Li(Si)), V e L^{[to,tiIL^{S')), 
M>0, |V^|<1, MSi,MS2 eLi([To,Ti]xS^). 
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2. The Einstein equations i3.27i , i3.2Si hold in the distributional sense. 

3. For every convex weak entropy flux If" to the Euler equations, the entropy inequality 

v„j« < g 

holds in the distributional sense. 

4. The following initial data are assumed in the distributional sense: 

(A„ A)(to, ■) = {A, A), {Ur, U)(to, •) = {U,U), 
(r]„ r]){To, ■) = {Ti Ij), (R„ R)(to, •) = {R, R), 



M(to,-)=M, V{tq,-)^V. 

The notion of distributional solution is analogous to the one already used in the previous 
section; see in particular the statement l|4.141 for a specific statement. It is important to observe that 
our regularity conditions are sufficient (and essentially necessary) to define in the distributional 
sense, all of the terms appearing in the Einstein-Euler equations. 

Monotonicity property of the area function 

From now on we consider a finite energy solution in the sense of Def inition IS . 1 l and we investigate 
its geometric properties. In particular, by definition, the area function R remains bounded and 
boimded away from zero. We consider the largest region D'^^^^^ of the future Cauchy development 
of the given initial data set which can be covered by a single chart in conformal coordinates. 

Specifically, following earlier work on classical solutions, we investigate the properties of the 
gradient VR and, by combining the evolution and constraint equations satisfied by R, we establish 
that second-order derivatives of the area function R belong to L!^(Lg), which, in turn, implies that 
the gradient VR is continuous. 

Proposition 5.2. The area function R satisfies the following properties: 

1. The functions Rr, Re are continuous in both variables and satisfy one of the following alternatives: 

Case (a) : Rr < -|Re| for all x, 6, 
Case (b) : Rr > \Re\ for all z, 6, 

Case (c) : R is constant and the spacetime is flat and vacuum. 

2. In Cases (a) and (b) the gradient of the function R is always timelike. 

3. In Case (a) the components of the gradient VR are uniformly controlled in L°°: 

IRel < -Rx <2sup|R,|. 

SI 

4. The second-order derivatives Ret and Rgg are uniformly controlled in L'^L^: 

sup r (|R,e| + |R00|)(T,-)|d0< r (|R.| + |Rel)(To,-)|d0+ f p{To,-)de. (5.1) 

T>To JSI Jfil 

In Case (a) (Case (b), respectively) the spacetime is contracting (expanding, resp.) in the future 
timelike direction. Our primary interest in this section is on the case of contracting spacetimes; 
expanding spacetimes were handled directly in areal coordinates, as discussed in the previous 
section. 
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Proof. 1. We generalize arguments due to Chrusciel |7l in the case of smooth vacuum spacetimes. 
Consider the functions R'^ := R-: +Re- We obtain 

R'g^ + {rjr-ile)R'^=F^ 

:= R{Ur ± Uef + 4^(^t ± Agf + _Re2(''-")(p + f). 

These equations, by assumption, hold in the sense of distributions. However, let us examine 
the regularity of the various terms involved, as follows. Our regularity assumptions show that 
(It - Tie) and F- belong to L'^{Lg), so that from our assumption that K'=^ e L;^(Lg ) we deduce the 
stronger regularity property 

R'^ e C(W^'i). 

In particular, both R"^ and, therefore, Rt,Rg are continuous in 6. This observation justifies the 
following calculations. 

Fix some time z. Take any point 6i e at which the function R{t, •) achieves a local extremum 
value; then, since the function Rg is continuous, we must have Rg{T, di) = and, therefore, 

(K'-R'+)(t,0i)=R,(t,0i)2>O. 

The case of equality is studied as follows. 

For almost every time t and every do e [0, 2n], we can write 

dd ^ ^ 
where we have observed that, since |/| < p, the functions F- are non-negative. It follows that 

R'^(T, do) < R'^iT, do + In) /„"""(')-')«)(-'9')''e' 
=: R'^(t,0o)H(t), 

Recall that the functions have just been proven to be continuous and observe that H is a positive 
quantity. Therefore, from the two inequalities 

R'-(t,0o) < R'-(t,0o)H(t), R'^(t,0o) < R'^(t,0o)H(t), 
valid for all do, we deduce 

R'-(t,0o)K'^(t,0o)(1-H(t))2>O. 

We recover the fact that the product (R'"R'"^)(t, ■) remains non-negative in the interval [0, 2n]. 
However, if now R'''"(t, Oq) - (for instance) for some do, then from 

R'^(t, 0) eXc-'"''*^'^'''*'^') = f F+(t, 6') dd' 

and the fact that F"'"(t, •) is non-negative, we conclude that R'"'"(t, •) changes sign only once, going 
from non-positive to non-negative values across do, which is impossible for a periodic and 
continuous function, except in the special case that R'"'"(t, •) is constant in d and, in addition, 
F-(t, ■) vanishes identically. 

The above arguments show that the function (R'~R'''")(t, ■) is either positive on [0, 27i] or else 
vanishes identically. 

2. The functions R'"*" (and therefore Rt,R9) are continuous in both variables. Indeed, we already 
know that Rgg e L!^(Lg) and, in view of the equation satisfied by the function R, 

R„ = Ree+M^(R)eC(4). 
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This leads us to distinguish between three different cases for all times: 

Case (a) : R" < 0, R'^ < 0, 
Case (b) : R" > 0, R'^ > 0, 
Case (c) : R" = or R'^ = 0. 

Suppose that, for instance, R'~{t,-) = and thus R'"'"(t, •) - 2Rg{T,-). From the first part of 
this proof, we know that the function R'"'"(t, •) is either positive, or negative, or identically zero. 
But, since Rg{T, ■) is 27T-periodic and continuous, it must vanish identically. Consequently, both 
K'='=(t, ■) = and the function K(t, ■) is constant in 6. We also deduce from the expression of that 
this is the case if and only if Uj, Ug,Aj,Ag, p vanish at the time tq. Moreover, it is not difficult 
to check that the fimctions U, A must be constant and p be identically zero everywhere in their 
domain of definition in the future of the hypersurface t = tq. The only non- trivial component the 
is the function i] which then satisfies a linear wave equation 

rirr - rjee = 0, 

whose solution takes the form 7](t, 6) = /-{t + 6) + f+ (t - 6) with /± e W^'^ . By setting u±:= % + d, 
the metric takes the form 

g = e-^^ f_{u.)f4u+)du.du+ + e^^{dx + Adyf + e'^^R^dy^, 

where the coefficients A, U, R are constants. Hence, the metric is flat. 

3. Finally, excluding Case (c) of the proposition, we can compute the norm of the gradient of the 
fimction R and we find 

From the conditions in Cases (a) and (b) we deduce that the component Rt can not vanish and 
therefore by continuity must keep a constant sign through the spacetime. 
Moreover, in Case (a) we can integrate the equation in R 

{d, + de)R'- =M^>0. 

This yields the sup norm control of both R'-, and in turn a control of both R^ and Rg. 

4. The estimate \5.1\ follows directly from ll3.28l l, by noticing that the right-hand side of (|3.281 is 
uniformly bounded in L^(S^). The latter indeed is easily checked from the fact that Ug,AT^,Ag 
are uniformly bounded in L^(S^) while p is uniformly bounded in L^(S^). 

□ 

In Case (a), that is, when the spacetime is contracting in the future direction, it will follow 
from our compactness framework in the next section that solutions in conformal coordinates do 
exist up to T ^ oo, as long as the function R does not vanish. As in the case of smooth spacetimes, 
one can distinguish between two situations: 

(al): {U,A,i],R,M,V) is a finite energy solution to the Einstein-Euler equations in conformal 
coordinates whose domain of definition D^^^^ coincides with [tq, X S^. In this case, the 
conformal coordinates do cover the whole of the Cauchy development. 

(a2): {U,A,rj,R,M,V) is a finite energy solution to the Einstein-Euler equations defined on a 
region strictly contained in [tq, X S^. Here, R admits a constant limit value Rmax which 
determines the boundary of the spacetime. 
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6 Compactness properties 



In the previous sections, we have successfully determined the natural weak regularity assump- 
tions required on the initial data and enjoyed by the solutions to the Einstein-Euler system. The- 
orem |23] now follows from the compactness framework developed in [20J. The main property, 
stated below, is that any sequence of spacetimes satisfying the natural bounds on the geometry 
and fluid variables must subconverge in a suitably strong sense so that the limit is a space- 
time with the expected regularity. The compactness property holds for exact solutions as well 
as approximate solutions (which should be dealt with by allowing for suitably controlled error 
terms). 

Proposition 6.1 (Pre-compactness property for the Einstein equations under Gowdy symmetry). 

Consider a sequence of spacetimes with finite energy determined by the geometric coefficients 17*"^, A^''\ Tf''\ 
R^"^ in conformal coordinates (or U'-'''', A^"^ , v^"', a^"' in areal coordinates) together with the (rescaled) fluid 
components M^"^, V'"'. Suppose that the regularity conditions stated in Definition \5.1\ (or Definition \4.5i 
hold, with the corresponding functional norms of the initial data being uniformly bounded in the parameter 
a = 1,1, . . .. Then, the following pre-compactness property holds: there exists a subsequence (still labelled 
with the upper script a) converging in the natural functional spaces and the limit is a solution with finite 
energy to the Einstein equations. 

To conclude, we consider briefly the case of the Euler equations in the flat geometry, that is, 
in the simplified situation that the coupling between the matter and the geometry is neglected. 
Hence, we focus first on the fluid evolution governed by the Euler equations: 



where /,< > and V e (—1, 1) are the unknown fields, and k e (0, 1) is a constant. Equivalently, in 
terms of the rescaled density M := p/{l - V-^) we have 

(M{l + k^V^)) +(M{l + k^)v) =0, 

' ^ (6.1) 

(m (1 + fc2) V)^ + (M{k^ + V^fj^ = 0. 

According to our definition, a pair of functions (M, V) defined on some set Q c (fi, t2) x and 
satisfying M e J^'t^i^g) and V e L°° with V e [-1,1], is called an entropy solution to the Euler 
equations l|6.Hl if the entropy inequalities 

{9°)t + < 

hold in the distributional sense for convex weak entropy pairs (3*, ?^). Note that, in the proposed 
formulation, the velocity may reach the limiting values ±1. 
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